Inhomogeneous States in a Small Magnetic Disk with Single-Ion Surface Anisotropy 



V. E. Kiree\E and B. A. Ivano\Q 

Institute of Magnetism NAS of Ukraine, 36-B Vernadskii avenue, 03142 Kiev, Ukraine 

(Dated: February 1, 2008) 

We investigate analytically and numerically the ground and metastable states for easy-plane 
Heisenberg magnets with single-ion surface anisotropy and disk geometry. The configurations with 
two half- vortices at the opposite points of the border are shown to be preferable for strong anisotropy. 
We propose a simple analytical description of the spin configurations for all values of a surface 
anisotropy. The effects of lattice pinning leads to appearance of a set of metastable configurations. 
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The progress of nanotechnology permits creation of 
ensembles of fine magnetic particles (magnetic dots) of 
nanometer scale, see for reviewi Magnetic dots in the 
form of cylinders or prisms have been made of soft mag- 
netic materials like Co and permalloySi 3 ^^ or highly 
anisotropic materials like Dy and FePt, seeAS Magnetic 
dots and their arrays are of interest both in the basic and 
applied magnetism with potential applications including 
high-density magnetic storage medial 

Usually a small magnetic particle is considered as be- 
ing in the monodomain state with a homogeneous sat- 
urated magnetization (or Neel vector for antiferromag- 
nets). During the last few years it had been established 
that the distribution of magnetization within the dots 
made of soft magnetic materials can be quite nontriv- 
ial; namely, various inhomogeneous states resulting from 
the magnetic dipole-dipole interaction appear. In re- 
cent years interest in such states for submicron particles 
has risen significantly. A small enough non-ellipsoidal 
dot exhibits a single-domain nearly uniform magnetiza- 
tion state, either so called flower and leaf states AiSiiiii^ 
When increasing the size of the dot above a critical value, 
vortex state occursAiii^i^iSiii The main property of 
such states is the non-saturating value of the total dot 
magnetization, nearly zero for vortex states and nearly 
saturated, but smaller that saturated, for leaf and flower 
states. 

O'Shea and coworkersi&iiLSI have observed non- 
saturated states for the rare-earth ferromagnetic gran- 
ules with high anisotropy and the size about of 5nm. 
A possible explanation of this fact is that these particles 
are in non-uniform states^ On the other hand, it is clear 
that the concepts of non-uniform states referred to above 
and caused by a magnetic dipole interaction cannot be 
applied directly to such small particles made with highly- 
anisotropic material. In this concern, some other sources 
of non-uniformity need to be found. 

The appearance of non-uniform states for small bee 
atomic clusters with taking into account the single-ion 
surface anisotropy have been shown numerically by Dim- 
itrov and WysiniSLS Garanin and Kachkachi in the re- 
sent workS investigated the effective anisotropy caused 
by such a non-uniform spin distribution for small mag- 
netic particles. The difference of the properties of the 
spins on the surface and in bulk could be considered as 



a defect destroying the homogeneity of a sample. It is 
clear that due to the surface a homogeneous ordering is 
distorted or even broken. 

In real magnets the surface could produce the surface 
anisotropy for two reasons. First, the main origin of mag- 
netic anisotropy can be caused by the anisotropy of spin- 
spin interactions (the case of exchange anisotropy). For 
this case even on an ideal atomically smooth surface the 
spins have different coordination numbers than in bulk, 
and consequently the intensity of the exchange interac- 
tion changes. For the surface exchange anisotropy the 
direction of the chosen axis is the same as in bulk and 
has no connection to the surface. This effect could lead to 
the non-uniform states in some special cases only, mostly 
in the presence of an external magnetic field, for example 
the surface spin-flop transitionf2i2^ and the states caused 
by the magnetic field for easy-axial ferromagnetsi2& Sec- 
ond, in real magnets surface atoms have a different en- 
vironmental symmetry. Thus, the surface distorts a 
crystalline field that acts on a magnetic ion, and the 
anisotropy is changed drastically. It leads to a specific 
single-ion surface anisotropy for the spins with a pre- 
ferred axis coinciding with the normal to the surface. 
This model is considered by Dimitrov and Wysin for fee 
iron clusters we would like to investigate this case 
both analytically and numerically. Note that the sur- 
face effects, in particular, the surface anisotropy, have 
been considered by many authors j2L2£ but in most of 
these works the ground states has been assumed to be 
homogeneous, and the surface terms are only accounted 
in dynamics. On the other hand, it is obvious that for 
fine magnetic particles the role of the surface becomes 
much more important than for bulk materials. The ef- 
fects cased by the surface considered as a defect are pro- 
portional to A" 1 / 3 , and their role increases when the size 
of the particle tends to the nanometer scales. 

Note that similar problems arise in the other domains 
of condensed matter physics, where a role of surface is 
important. These are textures in liquid crystals^ and in 
a superfluid 3 He, see^ For the A-phase of 3 He ( 3 He-A) 
the unit vector order parameter I, I 2 = 1 is perpendic- 
ular to the surface of a vessel. 3 He can not be in equi- 
librium with its own vapor; it fills the vessel completely 
at temperatures when it is superfluid (T < 2mK). Thus, 
the vector I should be perpendicular to the surface of the 
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3 He-A sample. The analysis shows that the order param- 
eter becomes non-uniform, and, moreover, it is singular 
for any simply-connected vessel^ 

It is clear that such effects may be observed in all 
finite samples of ordered media with vector order pa- 
rameter and a strong surface anisotropy of the form 
B(m ■ n) 2 , where n is the normal to the surface, m is 
the order parameter, and B is the constant of single-ion 
surface anisotropy, which orients m with respect to the 
surface. For the 3 Hc-A, the boundary condition could 
be described as a limit of an infinitely strong surface 
anisotropy B < 0, \B\ — > oo, with easy axis perpen- 
dicular to the surface. The concept developed for 3 He-A 
could be a good guide for a theory of fine magnetic par- 
ticles with surface anisotropy. On the other hand, the 
situation for magnets is more general: the magnitude of 
the surface anisotropy for magnets is finite, and the mag- 
netic moment could be inclined with respect to the axis 
of surface anisotropy. As we will show below finiteness 
of anisotropy could lead to the states with non-uniform 
spin distributions but without singularities. 

The outline of the paper is as follows. In Sec. [i] 
we discuss classical models for a small magnetic parti- 
cle supporting simplest non-uniform spin distribution, 
caused by surface anisotropy, which is planar and two- 
dimensional {2D) model. This means that the spins par- 
allel to one plane and the spin distribution depends ef- 
fectively on only two space coordinates (x,y). Sec. ITT1 is 
devoted to the planar continuum 2D model in the limit 
case of the infinite surface anisotropy, where exact so- 
lutions are found and analyzed. In Sec. II I II the same 
model will be considered for the case of finite anisotropy. 
Sec. IIVI contains results of direct numerical simulations 
for the 2D lattice models and the consideration of pin- 
ning effects that can be estimated from the continuum 
model. The analysis of thermal and topological stability 
is also done in this section. The last Sec. contains the 
resume of obtained results and a short discussion them 
in concern with other similar systems. 

I. MODEL 

There are two approaches to the analysis of the static 
and dynamic properties of magnetic materials: discrete 
microscopic and macroscopic. The microscopic approach 
is based on a discrete spin Hamiltonian in which the spins 
Si (quantum or treated quasi-classically, as will be done 
below) are specified at the lattice sites i. In discrete 
models the magnetic anisotropy can be introduced in two 
different ways: as single-ion anisotropy, and as anisotropy 
of the exchange interaction. To describe them, the spin 
Hamiltonian is chosen in the form 

h = JaSfsf + £ K 4 S ") 2 + E B °f><t)s? sf . 

<ij> i i 

(i) 

Here Sf is the projection of a classical spin on the sym- 
metry axis a of the bulk crystal. The summation in the 



first term is over all the nearest neighbors in the lattice, 
J a is an anisotropic exchange tensor. The constant K a 
and function B a p(i) describe the volume and the sur- 
face single-ion anisotropy energies, respectively. For the 
crystals with rhombic or higher symmetry, all tensors 
describing volume characteristics can be diagonalized si- 
multaneously. The tensor function B a p{i) is nonzero only 
near the surface and abruptly decreases in the depth of 
the sample. The surface creates another chosen direction, 
a normal to it, and enters a local system of coordinates, 
in which the tensor B a p(i) is diagonal. We neglect in the 
Hamiltonian Q a dipolc-dipole coupling and a Zeeman 
interaction with an external magnetic field. 

We shall use a simple version of Q with an uniaxial 
symmetry for the bulk properties (2 as a chosen axis, for 
definiteness) and with nearest neighbors interaction only: 

i,S 

+ K z J2(S;f + Bj2(S,n) 2 . (2) 

i i' 

Here J is the exchange integral, A is the anisotropy pa- 
rameter of the exchange interaction, and S are the vec- 
tors of the nearest neighbors, the summation over i' in 
the last term includes only the surface sites, where the 
number of the nearest neighbors differs from the volume 
one. To more adequately compare the lattice and con- 
tinuum models, we assume that the vector n is a normal 
to the surface, but not a direction given by the Miller 
indices. 

The sign of the exchange integral plays no role for the 
statics of non-frustrated magnets with a bipartite lattice. 
Moreover, a model without dipole-dipolc coupling is more 
adequate for antiferromagnets than ferromagnets. For 
simplicity we use below the ferromagnetic representation 
of spin distributions, i.e. J > 0. The transition to the 
antifcrromagnetic case for a bipartite lattice is trivial: 
we introduce sublattices and change the directions of the 
spins in one of them. 

The continuum approximation of is based on a free 
energy functional W[m] that depends on the local nor- 
malized magnetization m(r), m 2 = 1. Using the stan- 
dard smoothing procedure of a lattice model, we write 
down the functional W[m] as 

W[™] = ~ / ^f{ja 2 [(Vm x f + (Vm y f 
L Jn a 

+ A(Vto 2 ) 2 ] + Km 2 + aB{mn) 2 S(r - r s )\ . (3) 

Here f2 is the volume of the particle, the vector r s pa- 
rameterizes the surface, S(r) is the Dirac delta- function, 
a is the lattice spacing, and S is the cross-section area. 
The solution of the Euler-Lagrange equation for © gives 
spin configurations with a preferential direction close to 
the surface. It is clear that the measure of inhomogene- 
ity depends on the problem parameters and the sample 
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shape. A simple consideration shows that for the fixed 
shape there are only two relevant parameters. The first 
one is the characteristic radius R/tq, where tq = ay J/K 
at A = 1 or r = 2ayA/(l — A) at K = is the magnetic 
length, defined in the same way as for bulk materials, 3 - 
The second parameter is the ratio of the exchange inte- 
gral to the surface anisotropy B/J. To simplify analysis 
we consider a model with a purely planar spin distri- 
bution. Such distributions appear for magnetic vortices 
at strong enough easy-axis anisotropy, A < A c , where 
A c ~ 0.7 when tq ~ ami In this case we obtain the one- 
parameter model characterized by the ratio B/J. As 
we will see below such a model demonstrates a wide set 
of inhomogeneous states and allows complete analytical 
and numerical investigations. We restrict ourselves the 
case of one more simplification, namely, a model which 
allows 2D spatial spin distributions, i.e. such distribu- 
tions which depend only on two spatial variables, say x 
and y. Apparently such a simplification is applicable to 
an island of a magnetic monoatomic layer shaped as a 
disk. For numerical simulation we will choose a frag- 
ment of the two-dimensional square lattice in the form 
of a disk. However, applicability of obtained results is 
not limited by this concrete case. It is easy to imagine 
situations when the same spatial-two-dimensional distri- 
bution is realized. As an example one can regard a ferro- 
magnetic particle with the volume easy-plane anisotropy, 
having a form of a cylinder with the base parallel to the 
easy-plane (the xy-plane) and with the axis along the 
z-axis. If one considers that the surface anisotropy con- 
stant B in is positive then the normal to the surface 
is the hard axis of the surface anisotropy. It is clear, that 
any planar spin distribution with S z (i) = ensures both 
the minimum of the volume and the surface anisotropy 
on the upper and bottom cylinder surfaces. In this case 
non-uniformity is caused only by the lateral cylinder sur- 
face, and one can expect that the distribution will be a 
spatial-two-dimensional one, with the same character as 
for the purely two-dimensional problem. 



II. A STRONG BORDER ANISOTROPY IN A 
CONTINUUM APPROACH 

We shall start from the simplest model to describe 
effects of surface anisotropy. Consider a disk-shaped 
(or cylinder-shaped, see above) magnet, with ccy-plane 
as an easy-plane, and assume a 2D spin distribution. 
We assume that the magnetization is a two-dimensional 
unit vector, in a polar mapping: m z = and m± = 
e x cos 4> + e y sin cj>, where (e a , e y , e z ) is the basis in the 
spin space and <f> — <p(x,y) is the angle between m and 
e x . The magnetic energy of the disk takes the form: 



W[</>] = JS 2 



1/dSCv, 



b j d X cos 2 {<p-x) 



(4) 



are the polar coordinates in the plane of magnet. The 
parameter b is proportional to the constant of a border 
anisotropy, b — (B/J)(R/a). We choose b > 0, and the 
preferential surface directions are tangent. This choice 
is motivated above; one more reason is that such an ef- 
fective term can be used to model the magnetic dipole 
interaction; 3 — The function 4>{p, x) may have singularities 
inside the disk Q. Minimal configurations for the energy 
(0} are constructed from solution of the respective Euler- 
Lagrange equations, which is the scalar Laplace equation 



V 2 <t> = , 
with the boundary condition at p = R 



op 



-bsm2[(f>(p,x) - x] = 



(5) 



(G) 



Thus, this is a problem with a nonlinear boundary con- 
dition. 

In the absence of the boundary anisotropy, 6 = 0, ho- 
mogeneous solutions 4> = const satisfy simultaneously 
and ©, and this trivial case is not considered. First of 
all, we analyze possible solutions in the limit of strong 
border anisotropy, b — oo, when the problem becomes 
linear and can be solved exactly. The boundary condi- 
tion leads to the two possible solutions 4>(R, x) — X = ' =7r /2- 
Such ambiguity of the boundary conditions here differs 
from the classic internal Neumann problem of mathemat- 
ical physics and the relevant physics will be discussed 
below. The solutions in both cases can be constructed 
via harmonic functions, as well it can be done in two- 
dimensional electrostatics^ The general solution of the 
Laplace equation (j) can be written via a complex poten- 
tial u(z) of integer charges qk placed at the points Zk- 



(f> = Im[u(z)], u(z) ~ qk ln(z — Zk) + const . (7) 

k 

These charges have a simple physical meaning, they de- 
scribe well-known in-plane vortices, which have been re- 
peatedly discussed in regard to 2D magnetism. We intro- 
duce a complex representation for the coordinate plane 
xy, z = x + iy. The functional W[m] is rewritten as 



W[u] = JS 2 



where 





du 


2 


\ [ Pdp [ — 

2 Jo J\z\=p %z 


dz 




+ lf 


dz (e ^ i) 2 


8 J\z\=R 


iz 


e 



■ exp(u — u*) 



(8) 



(9) 



Here f2 is the area of our disk-shaped magnet with the 
radius R, the contour T is the border circle, and (p, x) 



In the continuum approximation the energy W[u] is 
logarithmically divergent close to points where the in- 
plane vortices (charges) qk are placed. To describe these 
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singular solutions in the continuum model we have to 
introduce a cut-off parameter of the order of the lat- 
tice spacing. Singularities cost much energy, and one 
could expect that configurations with a global minimum 
of W[u] should be sought among the nonsingular func- 
tions u(z) in the area f2 or functions with a small number 
of singularities. 



A. Vortex-like configurations 

The simplest s olutio n with one singularity is a centered 
generated by the functions u = 



vortex, see Fig. 1(a) 



\nz ± i-k/2 with the energy 



E„ = JS' z n\n{ — 



u(i) = ln(,-;--,,) + hi( z-??-\±i~ 



(10) 



where r e is a cut-off parameter for vortex states of the 
order of the lattice spacing a. Besides these solutions 
the others are non-centered vortices for infinite 6, see 
Fig. |l(b)| generated by 

where \zq\ < R , 
(11) 

with the vortex placed at the point zq. They also satisfy 
the conditions (j) = x ± 7r /2 on the border T. As seen 
from Hll[) . the interaction between the vortex and the 
border, which may be considered as a consequence of 
the boundary condition © , is equivalent to the coupling 
between the vortex and the image vortex placed outside 
the disk at the inverse symmetric point respective to the 
border circle. The calculation of the energy covers only 
the area and the singularity of the reflected charge 
gives no effect. The energy of non-ccntcrcd vortex for 
infinite b (fixed boundary conditions = x±7r/2on the 
border) is given by 



JS 2 ir 



lu| * 



In 1 



i*o r 

R 2 



(12) 



The first term coincides to a proper energy of the vor- 
tex given by ljlTj|> . and the second term is the energy of 
the interaction between the vortex and the border; it is 
a repulsive one. Besides it another force acts on the vor- 
tex: the vortex has a tend to escape from a finite area in 
order to decrease |Vm|, and, thereby being attracted to 
the border. In the case of b = oo the repulsive force pre- 
vails, the vortex is stabilized at the furthest point from 
the border, and the second term in (fP^ is absent. 
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(a) Centered vortex preferable for a strong border 
anisotropy. 
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(b) Non-centered vortex preferable for a weak border 
anisotropy. 



FIG. 1: Numerically calculated vortex-like states for the dis- 
crete model with A = 0, K z = and R — 5a. 



B. Configurations with two half- vortices on the 
border 

The above considered vortex-like distributions are 
some of the simplest spin distributions, minimizing sur- 
face anisotropy, not only for the circle shape but for a 



border in the form of any simple contour. Indeed, going 
around a simple closed contour, the normal n to it turns 
to 360°. This means that the topological characteristic 
of the planar unit vector, so called vorticity^ q equals to 
±1 for the vector n. Obviously, those magnetic vortices 
having the vorticity q — 1 are quite probable candidates 
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to realize the energy minimum. Nevertheless, vortices 
with any q ^ admittedly possess singularities inside the 
sample. The analysis of such distributions where magne- 
tization has no singularities in the bulk is of interest. A 
simple analysis demonstrates that in this case, as well as 
for 3 He-A singularities should appear on the border. 

To explain this, consider the behavior of the vector 
field m = e x co$<f> + e y smcf> on the border circle T. The 
boundary condition requires that the vector m be par- 
allel to the border. It can be presented by two ways: 
m may be parallel or antiparallel to the tangent vector 
r = n x e z . Assuming that m is nonsingular inside O, 
the circle T can be divided into an even number of alter- 
nating regions: in half of them m has to rotate clockwise 
and in the others — counterclockwise. Thus, besides the 
above considered vortex-like solutions, there exist config- 
urations regular inside the circle Q and with singularities 
on the border, see Fig. 3(c) (Such singularities in the 



three-dimensional case are referred to^S as vortex lines.) 
The simplest two-singularity solutions can be written as 



l(z) = ln(z - Re^ 1 ) + ln(z - Re i<t>2 ) + iir±i{ 



(13) 



This is a field created by two charges placed at the bor- 
der points Re 1 ^ 1 and Re 1 ^ 2 . It is easy to check that the 
conditions = j;±it/2 are satisfied on the border wher- 
ever where <f)(R,x) is defined. To calculate the energy 
thoroughly we have to introduce the cut-off parameter r' e 
and integrate over the disk f2 except two half-circles of 
radius r' e centered at the charges. Under the condition 
that the cut-off regions do not overlap, R((f>i — 02 ) ^> a, 
the energy of the configurations are 



E\ lv = J S 7T 



lnl ^ ) -ln2-ln 



sin ■ 



(14) 



Here r' e is the corresponding cut-off parameter. The con- 
tinuum approximation does not provide a relation be- 
tween r e and r' e and we used numerical calculations for 
the lattice model to find it out. These calculations show 
with a good accuracy that r e = r' e , and we will assume 
that in the following. The minimum of 114|) is achieved 
for charges placed at the opposite points of the border, 
it is given by 



hi 



R 



In 2 



(15) 



Thus, the interaction of surface charges with each other is 
also repulsive. Comparing the expressions (jlOJI and l)15|l. 
we see that the energies for both configurations are loga- 
rithmically diverged and differing by the constant. Thus, 
the configuration with two half-vortices at the opposite 
points of the border is preferable to the single vortex for 
Ay-model. 



III. FINITE VALUES OF A SURFACE 
ANISOTROPY 

In this section we consider the case of a finite surface 
anisotropy. At b < oo the boundary condition © is non- 
linear. It is easy to see that the only centered vortex from 
all configurations with the vortex inside the sample is an 
exact solution for any finite values b. A non-centered vor- 
tex is not a solution of our problem at finite b < oo. Such 
states are absent in the continuum model, but they be- 
came metastable in the discrete model because of lattice 
pinning. The numerical calculations shows that their en- 
ergies depend weakly on the surface anisotropy constant. 
This class will be considered in Sec. II VI 

The solutions with two half- vortices on the border l|13|) 
for finite anisotropy b < oo transforms to non-singular 
solutions with two vortices placed outside the disk at 
the opposite points zq and — zo, where \zq\ > R. This 
distribution is generated by the function 



i(z) = ln(z — zq) + ln(z + zq) + in ± in/2 



(16) 



The particular exact solutions of the problem (JSJ), © 
with an arbitrary b have been found by Burylov and 
Raikher— for a distribution of the vector director near 
the surface of a cylindrical solid particle embedded in a 
monodomain nematic liquid crystal. Using of the bound- 
ary condition © for the function l|ltj|) gives the value of 
Zq in the form 



|z | 2 -i? 2 (l + V / T+^)/& 



(17) 



For such values zq the boundary condition 10 satisfy 
exactly. The energy of the configuration is equal to 



E hv {b) = JS 2 it 



In 



2JS 2 



b- 



1 + VT+52 



x tan a; dx . (18) 



The latter term arises due to the cut-off close to the half- 
vortices which are introduced for |zo| — R < r e , and 4>q = 
arccos[(|zo| — R)/ r e]- Its contribution is important for a 
high enough surface anisotropy, B > J only, see Fig. |3(c)| 
It is easy to see that for any finite b the energy of the 
two-charge configuration is lower than its limit value 11511 . 
and it decreases monotonically with decreasing b. An- 
other limit case of small surface anisotropy 6-^0 leads 
to the almost homogeneous distribution m, see Fig. 3(a) 
with the nearly zero energy E = JS 2 irb. When b in- 
creases, the vector field m is curved to the diametrically 
pair of points, and the energy increases, see Fig. |3(c)| 
These features are in good agreement with that obtained 
numerically for discrete finite system. The dependency 
of -Ehv versus the surface anisotropy B from the Hamil- 
tonian J2J) is plotted in Fig. [5] together with that for the 
continuum model 10}. The discrepancy of the curves is 
connected with the discreteness effects, which are impor- 
tant for small samples, for larger system radius (the value 
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B/J 

FIG. 2: The dependence of the energy of the minimal con- 
figuration versus the surface anisotropy for the disk with the 
radius R = 10a. The thin line is a two-charge approximation; 
the thick line is a numerically calculated result for the lattice 
model. 
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(a) Almost homogeneous configuration, B/J = 0.2. 



of R till R = 30a has been used). In the case of b — > oo 
the vortex energy is higher than the two-charge config- 
uration energy for XY-model, and the vortex states are 
also metastable for any finite b. 



IV. NUMERIC SIMULATION AND LATTICE 
EFFECTS 

For our model with rather strong volume and surface 



anisotropy, the characteristic size is | Vm| 



and it is 



not obvious that effects of discreteness can be neglected. 
An exact analysis of the discrete model requires numer- 
ical calculations, but some qualitative results can be ob- 
tained using the lattice potential method. For a direct 
numerical simulation we basically used the XY-model 
with A = 0, i.e. with an extremely strong easy-plane 
anisotropy (some results concerning the finite A will be 
discussed in conclusion). 



A. Numerical simulation 

For numerical calculation of the equilibrium states we 
started from the discrete Hamiltonian for the magnetic 
energy J2J). Calculations have been performed starting 
from a random initial configuration or from a configura- 
tion given by (JJJ with constants Zk, q appropriate for a 
considered problem. The energy minimization has been 
performed through a Seidel-like algorithm with the suc- 
cessive exact solution of the local equilibrium equation 
for a fixed site that can be obtained from the following 
one-site energy 



D 
~2 



E L =-SHy-(Sn) 2 -^S 2 



(19) 



where \i is a Lagrange multiplier for the condition \ S\ = 1 
and H = JJ2 K> (S x e x + S y e y + XS z e z ) is the effective 
field created by the nearest neighbors of the fixed site. 
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(b) The charges are far from the border, B/J = 0.5. 
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(c) The charges are close to the border, B/J = 2.0. 

FIG. 3: Minimal non-topological configurations for the dis- 
crete model JSJ with A = 0, K z — and R — 5a for different 
values of border anisotropy. 
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FIG. 4: Metastable configurations for the R — 8a disk. Shown 
approximately 3 • 10 4 dots. Regions "1" and "2" are shown in 
Fig. and El 



The term with B is present for the border spins only. In 
a simple case B — we obtain S = H/\H\. When B ^ 
a more complicated analysis of the roots of the equilib- 
rium condition dE^/dS = is needed. Among these 
roots S mm we choose the value that gives the deepest 
minimum of E^. For all minimizations we observe that 
this procedure converges to one of the stable configura- 
tions. The configuration appearing during minimization 
energy process was mainly dictated by the choice of the 
initial configuration. 

In order to explore all metastable minimal configu- 
rations in the lattice model J2J with A = we per- 
formed more than 10 minimization procedures accord- 
ing with the described scheme. Initial configurations and 
the surface anisotropy B are chosen randomly. The ob- 
tained energy values are presented by dots on the plane 
(B/J, E/JS 2 ), see Fig.0J the system size is chosen small 
enough to show a discrete nature of the possible states. 
Such an analysis allows to judge both the energy absolute 
minimum for a given B/J and the presence of metastable 
states. It is seen that in some plane regions (marked as 
"1" and "2") the dots are grouped in more or less well- 
defined lines, which obviously corresponds to the most 
stable states and describes the dependence of their en- 
ergies on b. The characteristic regions are present in 
Figs. 03 El The region (marked as "3" ) , in which the 
dots are distributed practically randomly (in fact, there 
the dots also are fitted by lines), corresponds to high 
energy states. They are not subjects of interest. To 
classify the spin states the positions of singularities of 
the function <fi(x,y) have been analyzed numerically and 
the positions of poles (vortices), which are placed inside 
of a disk or on its border, have been obtained. Such 
an analysis demonstrated the presence of all states de- 
scribed above, including non-centered vortices and states 
with non-symmetrically placed surface singularities, but 
yet some less favorable states namely antivortices with 
the distribution like cf> = —\ + const, where x is the po- 
lar coordinate, instead of 4> = \ + const, characteristic 
for vortices. Let us discuss the obtained results. 
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FIG. 5: Minima with two half-vortices for the R — 8a disk. 
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FIG. 6: Vortex minima close to B = for the R — 8a disk. 



First of all, the given analysis has confirmed that 
the symmetrical states with two singularities possess the 
minimal energy. In the region of the small anisotropy 
and energy E < 4.0 (here and after energy values are 
presented in units of JS 2 ) only state with symmetric half- 
vortices are present, see details of this region in Fig. [SI At 
B/J> 0.5 other well-defined lines of dots appear, which 
also correspond to states with two half-vortices, however 
with broken symmetry. These states have higher energy 
and they are unstable at small B, but at larger B they 
become metastable due to surface pinning effects. With 
B increase first of all the states pinned in the vicinity of 
non-regular regions of a surface, which result from cut- 
ting a circle specimen from the square lattice appear. 
With further increase of B/J > 1.5—2 the number of 
asymmetric states grows. 

The second interesting region of plane at the energy 
E ~ 10 corresponds to vortex states. Its details at small 
B/J are depicted in Fig. It is worth to note that ac- 
cording to the analytical consideration the centered vor- 
tex presents at any B and its energy does not depend 
on B. Besides that state there exist non-centered vor- 
tices stabilized by the lattice pinning. Since the state 
with non-centered vortex at the finite b is not an exact 
solution (unlike to the non-singular case), they will be an- 
alyzed numerically in the next subsection with a simple 
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FIG. 7: The energy of a non-centered vortex versus its dis- 
placement for three values of surface anisotropy: (a) B/J = 
0.4 > B c (attraction to the center), (b) B/J = 0.04 ~ B c 
(equilibrium), (c) B/J — 0.004 <C B c (repulsion from the 
center). Radius of the disk is 30a. The radius of the pinning 
region (R p ~ 24) is maximal for the largest value of B and 
decreases for lower values. 



qualitative model of pinning. At large enough anisotropy 
the energies of non-centered vortices are higher compared 
to the state with the centered vortex. However, at small 
anisotropy, (B/J < 0.2 for the system size R = 8a used 
for Figs. 010 an interesting effect emerges: non-centered 
vortices become more favorable than the centered one 
(the correspondent region marked as "a"). This effect 
could be described as a change of the sign of the effective 
interaction between the vortex and the border at some 
value B = B c . (Let us remind that the case B = oo 
corresponds to fixed boundary conditions, while the case 
B = corresponds to free boundary conditions, which 
are associated with repulsion and attraction of the vortex 
to the border, respectively, see^i.) This effect is present 
also for big values of the radius, see Fig. [7| in which is 
plotted the vortex energy calculated in the model J2J) ver- 
sus its displacement for three value of surface anisotropy 
and the radius R = 32a. The characteristic value of 
the surface anisotropy B c decreases inverse proportion- 
ally to the system size. The values found numerically for 
R = (5 -r 30)a can be extrapolated by the dependence 
S c /J~ \.2{a/R). 

At B — the vortex and antivortex have the same en- 
ergy and in the region of extremely small B antivortices 
are also reliably observable, see Fig. region "b" . How- 
ever, when B increases the energy of antivortices grows 
rapidly and we do not discuss them. 



and the pinning effects. The simplest way to describe an- 
alytically lattice effects and, in particular, to investigate 
the local stability of metastable states, is to introduce an 
effective periodical potential (Peierls-Nabarro potential) 
into the continuum model. Schnitzer shows^ see also^i 
that for in-plane vortices this potential is independent 
of the values of out-of-plane anisotropy parameters (for 
A < 0.8) and can be presented in the simplest form as 
Upn{x, y) = nJS 2 n[sm 2 (xTr/a) + sin 2 (j/7r/a)], where the 
origin is chosen at the point which is equidistant from 
lattice sites, and the numeric parameter k ~ 0.200. 36 
The potential minima are attained at all points like 
r = ne x + me v , where to, n are integers, \e. x \ = \e y \ = a, 
and the saddle points are at (n + l/2)e x + me y and 
ne x + (to + l/2)e y . A metastable state with a vortex 
shifted from the center to the point r exists only when 
the sum E(r) — E v (r, b) + UpN(r) has a minimum at this 
point. The loss of stability manifests itself as ruptures of 
lines in Fig. |3] and El see also Fig.[7| 

Then it is easy to show that the non-centered vortices 
are held by the pinning potential and are stable if their 
coordinates are inside the circle of radius R p . The radius 
of the pinning region R p is determined from the explicit 
expression (|f 2fl for the energy of the vortex placed at the 
point ro as 



dE (d) 



(ro) 



dr 



r =R P 



UK 

a 



(20) 



and the case a< J? leads to R p = R — aj k. 

Thus the vortices can be pinned everywhere inside the 
sample except the thin strip close to the border. Their 
energies relative to the zero level of the centered vortex 
lie in the band of the width ~ J\n(R/a). Such states 
are frequently observed in numeric simulations for the 
discrete model when initial configurations for the mini- 
mization are chosen randomly. 

Although a detailed analysis of thermal fluctuations 
and decay of metastable states is beyond the scope of 
this work, their role can be discussed on the basis of 
the previous estimates. The above introduced R p is the 
radius of the region where pinning disappears, i.e. at 
r — > R p the barrier height separating states with a vortex 
placed in adjacent lattice sites, becomes to zero. It is also 
reasonable to introduce the function R p (E), such that at 
r < Rp(E) the barrier height between these two states 
is higher than some value E. Naturally, R P (E) 



77inax 



where iJ™ ax 



E -► 0, Rp(E) -► at E 
is the maximal pinning energy. For intermediate region 
E -C E™ ax a simple calculation yields 



Lattice effects and vortex stability 



R — R p (E) 



a £ b max 

- £max _ E > 



(21) 



For non-uniform states the lattice pinning of singular 
points in the spin distribution both vortices and surface 
singularities (half-vortices) plays an essential role. It is 
interesting to discuss such points in more details. The 
continuum model neglects a discrete nature of crystals 



and for all values of £?™ ax 



E - £™ ax the value of R P {E) 
is again near to R. Thus, the role of thermal fluctuations 
at kgT <C £™ ax can be considered as negligently small, 
and the above described metastable states may be man- 
ifest as long-lived ones even for finite temperature. On 
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the other hand at fc^T > £?™ ax metastable states like 
the non-centered vortex will not be manifest and only 
the centered vortex should be considered. 

For two-charge configurations the lattice potential also 
creates others metastable configurations with higher en- 
ergies than the energy of configurations with maximally 
separated charges. Their analysis is similar to the one 
that has performed for the case of a non-centered vortex. 
Two pinned charges on the border can be approached 
only down to the angle <p p = \<fii — 2 | — na/nR. Con- 
sideration of thermal fluctuations can be done for non- 
centered vortices as well and it leads to the similar re- 
sults, practically all such states are metastable. 

In conclusion of this section discuss the stability of 
vortices as a topologically nontrivial configuration under 
transform to non-topological one. Inside of two topologi- 
cally different classes of states — with vortex or with two 
surface singularities — effective relaxation to the most 
favorable state inside of the given class is possible. How- 
ever, the previous results show that the vortex-like con- 
figurations with the centered vortex have higher energy 
than the two-charge configuration, and are metastable. 
Therefore, the state with centered vortex may relax to- 
ward the most profitable state with two surface singu- 
larities. The simplest scenario of the vortex decay is the 
following. The vortex moves to the nearest point on the 
border and its counterpart moves also to it. The point, 
where they merge, is a saddle point of the path with the 
energy i5 sac i = JS 2 ir ln(i?/e) over the centered vortex en- 
ergy. This state is referred as a boojum or fountain in 
the 3 He theory) 3 ^* 3 ^ where it is a true minimum. Fur- 
ther, the merged charges decouple and move along the 
border: one — in the clockwise direction and another 
— in the counterclockwise direction to the most distant 
positions. Thus, the energy i? sa( j is the barrier height 
between the two classes of configurations, and it can be 
used for the analysis of a thermal (or quantum tunnel- 
ing, for low temperature) decay of vortex states. Note, 
the barrier is nothing to do with pinning potential. Its 
value does not contain the parameter k, but it is propor- 
tional to ln(i?/a) and is much higher that the exchange 
energy JS 2 . Thus, vortex states can be stable even at 
high enough temperature comparing with the Curie tem- 
perature T c ~ JS 2 , and the probability of the decay of 
the vortex state is very low even at the temperatures 
comparable with T c . 

V. CONCLUDING REMARKS 

A strong surface anisotropy for easy-plane Heisenberg 
magnets destructs the homogeneous ordering and leads to 
the two types of static structures: the vortex state and 
the state with pair of half-vortices on the surface. For 
finite anisotropy the latter state becomes non-singular. 
This state is energetically favorable for all finite values 



of surface anisotropy. The energy gap between it and 
the vortex state is of the order of the exchange energy, 
but the energy barrier is much higher that the exchange 
energy. The strong bulk anisotropy leads to well pro- 
nounced effects of lattice pinning, and large number of 
metastable states appears as well. 

It is interesting to compare these results to those which 
have been obtained for fine particles made with soft mag- 
netic materials such as permalloy magnetic dots, where 
the non-uniform states are caused by the magnetic-dipole 
interaction. The common point for these cases is not 
only the presence of the vortex state but also the pres- 
ence of non-topological non- uniform states, leaf or flower 
statesp^£^U£ The distinction consists in the fact that 
for soft magnetic particles there are non-singular vortices 
with the out-of-plane magnetization component while in 
our problem with the strong bulk anisotropy the only 
in-plane vortices with a singularity are presence. It is 
likely that in virtue of this for permalloy particles there 
is a very much pronounced transition from the vortex 
state to the non-topological one with the system size de- 
creasing, while in our problem the vortex state is always 
less favorable energetically. It is worth to note that our 
preliminary numerical data indicate the appearance of 
such a transition at a weak easy-plane anisotropy; an ex- 
tended discussion of this problem is beyond the scope of 
the present work. 

It is also interesting to note that the spin distribution 
in the non-singular state of our 2D problem resembles 
the distribution having axial symmetry and the plane of 
symmetry perpendicular to the axis obtained by Dim- 
itrov and Wysin2L2£ for 3-D particles where both the 
volume and surface anisotropics are presented. Recently, 
the stable three-dimensional analog of vortices, hedge- 
hog configuration has been discovered for a ball-shaped 
particle with strong normal border anisotropy by numeric 
calculations. 38 On the other hand, for the superfluid 3 He- 
A, which is defined in terms of our model by use of the 
infinitely strong surface anisotropy and isotropic volume 
properties, the true minimum constitutes less symmetric 
state (boojum, or fountain) with one surface singular- 
ity and without the symmetry plane^2*2£ In our case the 
"boojum- like" distribution appears only for non-stable 
saddle point, which separates the vortex and non-singular 
states. 
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